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Abstract. The La Cantera Formation (Aptian) represents an underfilled lake system developed in an extensional basin during the Cretaceous 
hothouse and it bears an abundant and diverse fossil record. Our goal was to investigate paleoclimate conditions of the unit from a multiproxy 
approach. We analyzed two proxies: a) the stable isotope (carbon and oxygen) composition of basal Actinopterygii and Neopterygii fish remains 
and from indeterminate plant remains; and b) the clay mineral assemblage composition; providing two independent lines of evidence to 
reconstruct paleoclimate conditions including paleotemperature calculations for lake waters. The values obtained for fish remains for 6"°C 
range between -8.4 and -1.3 %o (6"C = x -5.57%o; 6 2.25); and "80 values that vary between -5.7 and -3.6%o (6'8O = x -4.33%o; o + 0.84); in 
the case of plant remains 6'°C values range between -25.8 and -22.8%o (81°C = x -25%o -4.33). Paleotemperatures obtained from 6780 for the 
water column ranged between 23.33 to 35.80 °C. Plant 5"3C signature allowed obtaining a 5'°Ccg3 composition of -5 %o. This agrees with the 
interpretation of the La Cantera originating from sediment deposition in an underfilled lake basin and with the global climate context of the 
Cretaceous. The clay mineral assemblage is dominated by illite-smectite indicating enhanced hydrolysis and seasonal rainfall supporting 
temporal warmup of the lake. Our data supports aridity conditions during the La Cantera paleolake existence. These findings suggest that 
geochemical proxies of paleontological datasets are accurate for paleoclimate reconstructions for the deep time records of the La Cantera 
lacustrine system. 


Key words. La Cantera Formation. Lacustrine. Underfilled lake. Stable isotopes. Fossil fish. Fossil plant. 


Resumen. REGISTRO PALEONTOLOGICO DE UN PALEOLAGO CRETACICO DE GONDWANA COMO PROXY PARA RECONSTRUCCIONES 
PALEOCLIMATICAS. La Formacion La Cantera representa un sistema lacustre underfilled desarrollado en una cuenca extensional durante la casa 
calida del Cretacico y posee un registro fésil abundante y diverso. Investigamos las condiciones paleoclimaticas para la unidad desde un enfoque 
multiproxy. Analizamos dos proxies: a) la composici6n isotdpica del carbono y del oxigeno de peces Actinopterygii basales y Neopterygii y de 
restos de plantas; y b) la composicion de los ensambles de arcillas; aportando lineas de evidencia independientes a la reconstrucci6n 
paleoclimatica. Los valores obtenidos para peces de 6"°C varian entre -8.4 y -1.3%0 ("C = x -5.57%o; o + 2.25); y para 6180 varian entre -5.7 y 
-3.6%o (680 = x -4.33%0; o + 0.84); en el caso de las plantas varian entre -25.8 y -22.8%o (6'°C = x -25%o; o + 4.33). Las paleotemperaturas 
obtenidas a partir del 5'°0 para la columna de agua fueron de 23.33 a 35.80 °C. La firma isotépica del 5"3C de las plantas permitio obtener la 
composicion de 6'Ccgz de -5 %o. Esto concuerda con la interpretaci6n de que La Formaci6n Cantera pudo haberse originado a partir de la 
depositaci6n de sedimentos en un paleolago underfilled y con el contexto paleoclimatico global del Cretacico. Los ensambles de arcillas estan 
dominados por illita-esmectita indicando un indice de hidrdlisis elevado y estacionalidad en las precipitaciones lo que apoya el calentamiento 
temporal del cuerpo de agua. Estos datos indican condiciones de aridez durante la existencia del paleolago La Cantera y sugieren que los proxies 
geoquimicos de las bases de datos paleontoldgicas son precisos para la reconstrucci6n del paleoclima local en tiempo profundo de la Formaci6n 
La Cantera. 


Palabras clave. Formacion La Cantera. Lacustre. Lago Underfilled. lsotopos estables. Peces fdsiles. Plantas fdsiles. 


SEDIMENTARY records (biotic and abiotic) are useful tools with fidelity archives in the sedimentary record for terrestrial 
high fidelity for paleoenvironmental and paleoclimate paleoecosystems, since they condense local and regional 


reconstructions, being lacustrine deposits the highest data that can be used as proxies of the environmental and 
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climate conditions of the past (Gierlowski-Kordesch & Kelts, 
1994). Such datasets have been applied successfully for 
Quaternary system reconstructions and deep time systems 
as well (e.g., Abell et a/, 1982; Wang et al, 2002; Benavente 
et al., 2021, 2022). Among the data preserved in lacustrine 
records, there are clay assemblages of the fine-grained 
deposits of lake basins that can be applied for such studies 
as well as the geochemical signature of stable isotopes (e.g, 
Benavente et al, 2019, 2021; Salduondo et al, 2022). The 
stable isotope signature in a lacustrine system can provide 
information about external (precipitation) and internal 
processes (stratification) (Leng et a/, 2005). Stable isotopes 
can also explore the different resources availability for the 
organisms that inhabit the system and how they exploit 
them (Longinelli & Nuti, 1973; Kolodny et a/., 1983; Grocke, 
2002; West et al, 2006; Cerling et al, 2007), like the food 
source or the niche they inhabit in the body of water 
(shallow marginal versus more profundal) (Pereira et al, 
2021). 

La Cantera Formation (Sierra del Gigante, San Luis 
Province) is a unit that represents a lacustrine system. Their 
sediments were deposited during Aptian times in an exten- 
sional tectonic context (Flores, 1969; Criado Roque et al, 1981; 
Rivarola & Spalletti, 2006). It bears a rich fossil record with 
high fidelity preservation that includes fish, plant (micro and 
macro), and insects (Puebla et a/,, 2012; Arcucci et a/, 2015; 
Castillo Elias, 2016; Giordano et a/, 2016). Even more, there 
is a study of the fossil record of the unit to infer paleoclimate 
conditions (Pramparo et al, 2018). Nevertheless, single 
proxy analysis has limitations and linked to the fossil record 
such proxy presents strong biases like preservation of the 
remains and other factors controlling their distribution, like 
local climate conditions instead of regional ones (Benavente 
et al, 2022). Therefore, coupling different proxies such as 
the geochemical composition of fossils and clay assemblage 
composition can be more accurate for environmental infer- 
ences. The isotopic signature of fossils is highly conser- 
vative in the sedimentary record and allows identifying 
environmental patterns that affected organisms in the 
system they inhabited (West et a/, 2006; Tutken et a/, 2006; 
Rey et al, 2020; Pereira et a/., 2021). For instance, the 
chemical composition of ganoid scales is generally used for 


robust paleohydrology and paleoclimate reconstructions 
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(lacumin et al, 1996; Grimes et al, 2003; TUtken et al, 2006; 
Amiot et a/., 2010; Pucéat et a/, 2010; Sisma-Ventura et al, 
2019). 

The La Cantera lacustrine system has been interpreted 
as an underfilled lake basin type (Castillo-Elias, 2016) ac- 
cording to Bohacs et al. (2000). This implies that it was 
characterized by an evaporative facies association, which is 
coherent with, among other diagnostic features, significant 
gypsum deposits found throughout the unit associated with 
mud cracks (Castillo-Elias, 2016). However, development of 
underfilled lake basins is not only controlled by paleoclimate 
alone but by hydrological factors (Benavente et al, 2019). 
This provides an interesting setting to reconstruct paleohy- 
drological conditions of such an evaporative system. 

The main goals of this study are: a) to test new geo- 
chemical proxies that can be extracted from the fossil record 
such as the stable isotope composition (carbon and oxygen) 
of fish and plant remains; and b) to obtain the clay assem- 
blage composition in order to provide a multiproxy interpre- 
tation for paleoenvironmental and paleoclimate conditions 


of the La Cantera lacustrine system during the Aptian. 


Geological Setting 

The San Luis Basin is located in the Southwest margin 
of the Sierras Pampeanas geological province, limiting to the 
Northwest with the Sierra Grande de San Luis Province 
(Yrigoyen, 1975). To the west, it extends up to the Bolson 
del Desaguadero (Mendoza Province) and to the north it is 
identified up to the Southwest limit of the San Juan Province 
(Fig. 1.1). The infilling of the basin is entirely continental and 
comprehends Aptian to Neogene deposits including the 
La Cantera Formation, that reaches up to 40 m in thickness 
at the basin center (Fig. 1.2). This unit outcrops at the Sierra 
del Gigante (Northwest of the San Luis Province) (Fig. 1.3) 
and represents a lacustrine system. These deposits are 
characterized by dominant ripple-cross sandstone and 
dark finely-laminated mudstone (Fig. 1.4) with abundant 
fossil remains (Flores, 1969; Pramparo, 1994; Castillo Elias, 
2016). 

The fossil record of the La Cantera Formation includes 
crustaceans (ostracods and conchostracans), insects, 
actinopterygian fish, algae, and micro- and macroflora 


(briophyta, monilophytes, equisetaleans, gnetaleans, and 
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Figure 1. 1, Geological map of the San Luis Basin. Modified from Arcucci et a/. (2015). 2, Stratigraphy of the El Gigante Group (Aptian—Albian), 
modified from Arcucci et a/. (2015). Pz, Paleozoic. 3, Outcrop photograph showing the Los Riscos (LR), La Cantera (LC), and La Cruz (Lcz) 
formations at the Sierra del Gigante deposits. 4, Field photograph of the La Cantera Formation at its type locality, where it is mostly constituted 
of dark, well-laminated mudstones interbedded with fine-grained sandstone. 


angiosperms) (Petrulevicius et a/, 2010; Arcucci et al, 2015; 
Giordano et al, 2016). The palynology assemblage has al- 
lowed to link the deposits to the Late Aptian (Pramparo, 
2012). The micro- and macroflora assemblage of Ephedripites 
sp., Cycadopites sp., Monosulcites sp., and Classopollis sp. has 
allowed also a single proxy paleoclimate inference of aridity 
indicators (Puebla et a/, 2012), with briophyta indicating 
local humidity (Puebla, 2009). 


MATERIALS AND METHODS 

A high-resolution stratigraphic log of the middle section 
of the La Cantera Formation was measured at 1:100 scale 
at the type locality “La Cantera de Gutierrez” (Sierra del 
Gigante) (32° 59’ 23.8" S; 66° 52' 52.6” W). Six levels were 
sampled and hand samples were collected for fossil extrac- 
tion and eight mudstone levels were sampled for clay as- 


semblage analysis. 


Materials 

Fish. Two groups of fish remains were analyzed with low 
magnification binocular microscope (Leica M80) at the 
fossil repository at the Universidad Nacional de San Luis 
(UNSL): 1) basal Actinopterygii up to 10 cm long and 2) 
ganoid Neopterygii up to 5.5 cm average of total length. The 
specimens were well preserved and they are housed at the 
Museo Interactivo de Ciencias (MIC), UNSL collection (Figs. 
2.13). 

For this study, uncatalogued specimens were used, since 
most of them were destroyed during sampling. Samples 
were taken from teeth and bones (Tab. 1) of nine basal 
Actinopterygii specimens (field number: C1, C2, C3, C4, C5A, 
C5B, C7, C8, and SLC2018) and from bones, fin elements, 
scutes, and scales of nine ganoid Neopterygii specimens 
(MIC-V707; and field number: T1, T4, T5, T6, T7, T8, T9, and TB2). 
Macroflora. Five samples of indeterminate stem remains 


have been analyzed (Fig. 2.4). 
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Figure 2. Images of fossils selected for sampling. 1, Articulated basal Actinopterygii. Scale bar= 5 mm. 2, Maxilar bone and teeth from a basal 
Actinopterygii. Scale bar= 1 mm. 3, Patch of scales from a ganoid Neopterygii. Scale bar= 1 mm. 4, Plant remains. Scale bar= 1 cm. 


Methods 

Carbon and oxygen stable isotopes. Carbon and oxygen sta- 
ble isotopes were sampled from fish and plant remains. 
These samples were analyzed on a Thermo Fisher Scientific 
GasBench Il, equipped with a PAL autosampler, and coupled 
to a ConFlow IV interface and a MAT 253 mass spectrome- 
ter (Thermo Fisher Scientific), to determine their carbon and 
oxygen isotopic compositions (613C and 6'80) at the Stable 
Isotope Ratio Facility for Environmental Research, University 
of Utah, Salt Lake City, Utah, USA. 

The samples were weighed using a Sartorious micro- 
balance and loaded into 4.5 mL flat-bottomed borosilicate 
vials (Labco) and capped with Labco butyl rubber septa. 
Then, the vials were flushed for six min on a PAL autosam- 
pler with ultrahigh-purity grade helium (99.999% He) ata 
flow rate of 50 mL/min. During flushing, the vials were kept 
in a heated (50 °C) aluminum block. The samples were re- 
acted with 10 droplets of phosphoric acid (H3PO,) (kept 
at 50 °C) to produce CO, gas. The samples were allowed to 
equilibrate for 12 h before analysis. During analysis, the CO, 
gas produced in the vials was collected using a sampling 
loop of 100 pL and then transported to the mass spectrom- 
eter. Nine individual injections per analysis were made for 
each sample, and for reference materials, and their average 
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was taken as the final number. Three sets of internal refer- 
ence materials were used to calibrate the system and our 
samples. Carrara marble, LSVEC, and Marble-Std were 
used as primary reference materials, and Marble-Std was 
used as a secondary reference material to cross-check the 
final number. Internal reference materials were calibrated 
against international standards NBS-18 and NBS-19. 
Values are reported as per mil relative to the Vienna Peedee 
belemnite (VPDB) standard. Precision was better than 
0.05%o for 5180 and 0.02% for 61°C. 

Mineralogical analysis. Mineralogical composition of bulk 
and the oriented 2 um size fraction of the eight mudstone 
samples (M1/18-M8/18) was analyzed through X-ray dif- 
fraction (XRD) analysis at the Centro de Investigaciones 
Geoldgicas in La Plata City, Buenos Aires Province, Argentina, 
using a Panalytical X'PERT Pro diffractometer equipment 
following Moore & Reynolds (1997). The eight rock sam- 
ples were disaggregated and powdered in an agate mortar. 
Clay minerals (Tab. 2) were obtained by pipetting the frac- 
tion <2 um ina suspension with distilled water following the 
Stokes Law. The samples were prepared according to the 
glass slide method: a) natural, air-dried sample in the labo- 
ratory at environment temperature; b) glycolated, sample 
exposed to the vapors of an ethylene glycol solution for at 
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TABLE 2 — Bulk rock mineralogy of mudstones from the La Cantera Formation, Aptian, Sierra del Gigante, San Luis Province, Argentina 


Sample Ka ll /S sm Ch An Qz Fs Pl Ca Py 
M8/18 XX X X X X x Xx 

M7/18 x XX X Xx x Xx Xx Xx Xx XX 

M6/18 Xx XX % x Xx Xx x XX 

M5/18 Xx XX X x x X Xx x X XX 

M4/18 Xx XX Xx Xx X Xx x XX 

M3/18 Xx XX X x Xx XX Xx Xx XX Xx 
M2/18 X Xx X x Xx X Xx Xx XX 

M1/18 Xx XX X Xx Xx XX x Xx Xx 

Abbreviations: An, analcime; Ca, calcite; Ch, chlorite; Fs, feldspar; Ka, kaolinite; Il, illite; I/S, interstratified illite/smectite; Pl, plagioclase; Py, 


pyrite;Qz, quartz; Sm, smectite; x, presence; xx, high intensity XRD peak. 


least 24 h; and c) calcined, the sample taken at 550 °C for 
two hours. X-ray diffraction patterns were obtained using 
Cu K-alpha radiation (ka =1.5403 A) operated at 40 kV, 40 
mA, and scanning speed of 0.04°/s, between 4° to 37° 28, 
for the total rock samples; among 2° to 32° for the natural 
samples of the clay fraction; 2° to 27° for the glycolated 
samples of that same fraction, and 3° to 15° for the calcined 
samples. The identification of main mineral phases and clay 
minerals was based on diagnostic peaks according following 
Moore & Reynolds (1997). 


Diagenesis assessment 

In order to determine the diagenetic overprint of the 
isotopic signature, the correlation coefficient between 
8180-93 and §'®Opo,, was analyzed for the samples (lacumin, 
1996; Michener & Lajtha, 2008), using the software 
STATISTICA (StatSoft, TIBCO Software, Windows). 


Paleotemperature calculation 

It was done following the equation from Longinelli & 
Nuti (1973) and modified by Tutken et a/. (2006) for fish 
teeth: 


Following the criteria from Benavente et al. (2021, 
2022), the 5'8O of modern meteoric water was obtained 


from Waterisotopes.org (http://wateriso.utah.edu/wateriso- 
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topes/pages/data_access/oipc.html) (Bowen & Wilkinson, 
2002), inputting latitude and altitude. Paleolatitude was 
calculated for the basin during the Cretaceous from the 
Paleolatitude Calculator for Paleoclimate Studies 
(http://www.paleolatitude.org/) (van Hinsbergen et al, 
2015) but calculation showed no significant difference with 
modern latitude; therefore, modern coordinates were used 
(32° 49’ S; 66° 59’ W). Paleoaltitude was inferred as ap- 
proximately 500 masl according to the development of rift 
lakes in the modern East Africa Rift (Bessems et al, 2008; 
Bergner et al, 2009; Moernaut et a/, 2010). 

We applied the Kolmogorov-Smirnov statistical test to 
the distribution of the data using Excel (Microsoft, 
Windows). Data distribution of both variables (6180, 673C) 
was normal; therefore, the Pearson correlation coefficient 
was calculated using InfoStat 2020 software (Universidad 


Nacional de Cordoba, Windows, Cérdoba). 


8°Cco2 calculation 
It was done following the equation from Farquhar et a/. 
(1989) simplified in Grocke (2002), assuming linear isotopic 


fractionation: 
513C giz (%0) = 5'?Cpiant + 20.22 


We also applied the equation proposed by Arens & 
Jahren (2000) to test the better fit to the dataset: 


5°Cco2 = (6°Cptant + 18.67)/1. 10 
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RESULTS 
C and O stable isotopes from fossil remains 

Fish samples from the La Cantera Formation have 6%C 
values that range between -8.4 and -1.3%o (6°Ca =X -5.57; 
o + 2.25); and 6'8O values that vary between -5.7 and - 
3.6%o (520 = X -4.33%0; o + 0.84) (Fig. 3, Tab. 1). Average 
values of 6"7C and 6"80 for both fish groups are very similar 
(8° Cact = X -5.91%0; 6 + 2.21; 6 °Cueop = X -5.26%o; 6 + 2.38; 
5 Onct = X ~4-48%0; 6 + 0.7; 8 '°Oneop = X ~4.2%o; o + 0.96) 
(Fig. 3; Tab. 1). 
Interpretation. Since both groups present similar averages 
in their isotopic signature, biological effects (metabolic or 
physiologic) are discarded. There are no significant differ- 
ences either when comparing the different materials ana- 
lyzed (bones, teeth, scales) (Tab. 1). One interpretation is 


that there is no difference in the fractionation effect for dif- 
ferent tissues in the organisms analyzed or alternatively 
perhaps, this effect has been homogenized by secondary 
imprint (diageneses and/or weathering) of the signature 
(Tieszen et al, 1983). 


Diagenesis analysis 

There is a low negative correlation for fish values be- 
tween 6780¢93 y §'8Opo,, with a correlation coefficient of 
r = -0.28 (Fig. 4.1). The correlation analysis by the taxonomic 
group of fish shows a similar result, with a correlation co- 
efficient of r = -0.41 for basal Actinopterygii (Fig. 4.2) and 
r = -0.39 for ganoid Neopterygii (Fig. 4.3) and also for the 
different tissues sampled (Fig. 4.4). 
Interpretation. The correlation coefficients obtained for the 
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Figure 3. Cross-plot of 6180 and 6"°C data of the fossil record (fish) of the La Cantera Formation (Aptian), San Luis Basin in the San Luis Province, 


Argentina. 
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complete fish remains dataset, each group separately, and 
even for different tissue types suggest that the stable iso- 
tope signature has been diagenetically altered (Kolodny et al, 
1983; lacumin et al, 1996; Tutken et al, 2006). Therefore, 
the fish remains dataset wouldn't be completely reliable for 


paleoclimate reconstructions. 


Paleotemperature 

Paleotemperatures obtained using 5'°Opo, values from 
fish remains vary between 23.33 °C and 49.43 °C (Tab. 1). 
Interpretation. Within the dataset, we consider that the 
warm end member presents unreasonable water column 
values for fish to inhabit (Grande, 2010). Considering the 
diagenesis analysis that suggests the primary isotopic 


signature might have been altered, we consider that the 
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most likely primary data set is the one resulting in a pale- 
othermal range of 23.33 °C and 35.8 °C. This implies dis- 
carding five values from samples from basal Actinopterygii 
and ganoid Neopterygii, because they were most likely sub- 
ject to diagenetic reset. Since this effect is potentially ob- 
served in different tissues and different stratigraphic levels, 
no clear diagenetic pattern can be inferred. 

Plant samples have 67°C values that range between -25.8 
and -22.8%o (86°C, = -25y; o + 4.33) (Tab. 1). 
Interpretation. The 6'3C signature of plant remains is very 
similar to modern C3 plants that show an average of -25.4%o 
and vary between -22.5%o and -27.9%o (Craig, 1954; 
Grocke, 2002). Our data are consistent with 6'7C wood 
values (Grocke, 2002). This suggests a high-fidelity preser- 


vation of the C stable isotope signature in this data set. 


Scatterplot: 5'8Oc(%o) T vs. 5'SOp (%0) T (Casewise MD deletion) 
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Figure 4. Linear regression graph of &'®Oc93 (%o) and 8'8Opo,, (%o) values using different samples, all belonging to the La Cantera Formation 
(Aptian), San Luis Basin, San Luis Province, Argentina. 1, All fish remains. 2, Linear regression graph of 618O¢93 (%o) and 6®Opg,, (%o) values of 
basal Actinopterygii. 3, Linear regression graph of 870.93 (%o) and &'®Opo,, (%o) values of Ganoid Neopterygii. 4, Linear regression graph of 
8"°0¢93 (%o) and &"®Opq,, (%o) values according to the tissue sampled from fish remains. 
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However, our sample size is limited to further constrain 
paleoenvironmental variables such as atmospheric CO> 
(Arens & Jahren, 2000). 


5°Cco2 composition 

Calculations made applying Farquhar et a/. (1989) equa- 
tion for 6"Ccgz present values that range between -2.6 
and -6.1%o with an average of -5.0%o. Calculations made 
following Arens & Jahren (2000) show 8'%C-q> values that 
range between -9.18 and -5.81%o with an average of -8.27%o 
(Tab. 1). 
Interpretation. According to previous reported values of 
5Cco2 for the Cretaceous, the formula that best fits for 
the data of the La Cantera wood remains is the Farquhar et 
al. (1989) (-3.0%o 61?C; Grocke, 2002), meaning that the 
atmospheric carbon isotopic composition was on average 


approximately -5.0%o 873C. 


Mineralogy 

Bulk mineralogy. XRD analysis shows the presence of 
quartz, plagioclase, calcite, feldspar, analcime, and pyrite 
(Figs. 5, 6; Tab. 2). Calcite is the most abundant in the 
samples followed by quartz. The analcime is recorded 
throughout the section, as well as plagioclase and feldspar. 
Pyrite is only recorded in level 1 (Sample M3/18). 

Clay mineral assemblages. Minerals identified in the mud- 
stones of the La Cantera Formation are illite, kaolinite, 
chlorite, smectite, and mixed-layer illite/smectite. Each 
level analyzed presents a different combination of these 
minerals. Illite is recorded in each level tested, whereas 
chlorite is present in seven of eight levels tested, except for 
the level 2 (sample M4/18). Kaolinite is recorded in the lower 
section, until level 3 of the stratigraphical section sampled 
(Fig. 5). Mixed-layer illite/smectite is recorded throughout 
the section in six levels (Figs. 5, 7) while smectite is recorded 
only in levels 3 and 9 (Fig. 5). 

Interpretation. All levels sampled contain illite. This clay 
mineral is the result of lixiviation and erosion, indicating high 
enhanced hydrolysis index due to high physical weathering 
(and low chemical weathering) compared to smectite or 
kaolinite (Chamley, 1989; Hillier et a/, 1995). The illite crys- 
tallinity is low, ruling out leaching and diagenesis processes 
(Fig. 7) (Fursich et a/, 2005; Barrenechea et a/., 2018). The 


presence of smectite and mixed-layer illite/smectite indi- 
cates that precipitation was seasonal (Singer, 1984; 
Chamley, 1989). 

There are no records of volcanic activity in the Late 
Aptian of the San Luis Basin, allowing to rule out volcanism 
as the origin for smectite (Do Campo et a/., 2010). However, 
we cannot discard the possibility of volcanic rocks in the 
source area as input for the smectite formation (Salduondo 
et al, 2022). The mixed-layer illite/smectite is interpreted 
as an indicator of leaching (Moore & Reynolds, 1997). The 
lower level, which records assemblages with illite and kaoli- 
nite, indicate enhanced leaching due to a high-water regime 
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Figure 5. Sedimentary log of the middle section of the La Cantera 
Formation (Aptian) with clay mineralogy and fish and plant samples 
in its stratigraphic location. I/S, mixed-layer illite/smectite; Sm, 
smectite; Kao, Kaolinite; Chl, chlorite. M, mudstone; Si, siltstone; S, 
sandstone; C, conglomerate. M1/18: clay samples. 
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plagioclase; Ca, calcite; A, analcime. 1, Sample M1. 2, Sample M2. 3, Sample M3. 4, Sample M4. 5, Sample M5. 6, Sample M6. 7, Sample M7. 


Figure 6. Total rock diffractograms of the La Cantera Formation (Aptian) mudstone samples in stratigraphical order showing: Q, quartz; PI, 
8, Sample M8. 


90 


BENAVENTE er. AL: GEOCHEMISTRY OF THE LA CANTERA FORMATION 


| ii m 
wo i | Clay M1-LC Seta 1) 00 Clay M2-LC Gheoated 
| I be | 
\ FW ; : 
A : ; 
2 | t a . bat © ea Wed A Ape a \ 
c Ny a i ea array s & | - tom fe 
s iy Yayhyrha fe 5 as gi 
ST ug, Mi . TF cael = ch bg Ad mdm | 
~ i wh | Mi naaenh a 1 . Wwe , A 3 ew Hp ra wry 
2 ; en a c = at wal 
this eH “yop ay Nh vif Anal e = ; . \. fe # | aa ‘ 
i A \) 8 [ 2898 | pSs ai 100 Wie eh | s 
ae a 88am St fe Ah | ge BB fice 
‘iayh | ie Ia yin V hee iby aye , \ hi t alien ff ih § J] 
AY \wy ie a me Woy th yee La Wy in we et Whe nintis 
iy is Pa N 
26 26 
i == Calcined \ —————_ icine 
<a } [clay Meic a) at | Wy i [clay Meio eyes 4) 
V'h Nay 
ts . Nw 2 
ay LA Sse : 
- bh, 8 m Ze 
2 \ 1 i i ‘wlan : F a z , MA a | : 3 \ 
k he “WA Fd ¢ I ; Wh? v hs = j ash h lt 
baal : 5 ih 1 ym 73 ary Na ai bea 88 | ‘ aha! weet 
at shu be Hn Ah, ol an ate oS i. I i ‘wy Pt : ‘ehh Nala Re 
Te TPT gg ye! ¢ 5 Wi ga 
ee ‘ 1 J Aya andy LN Ver ‘dah \ = cat I dada it M4 j } 
“out itty, pal aw My H Mat rien awe Aye eee) 
208 20 
P| i Clay M5-LC a) a ae Eon O) 
i a = tural 
\ Sa | BY 
> He i | 2 ; 
5 \ nt a i 3 ' 5 
= \ i Aufl Ng f 8 it = : & ag’ |\ 
hie Kh a e 8 wa fae J al iL i ' 
a OA sel 
a \y & ‘an Whi h i Vy Nw r sad | i wh nf ie 3 
"W i" hy! es Loe fk ‘ 8 "type PM at! ae | git 
ig Bes Ady Hg | eT, nail fie tah | 
y & é eect y ph +l Wa + ne ie wy 
VvaPlanc lf PO kant vl 
yey! allay mp 
26 26 
200 | ", Clay _M7-LC oyeoed| Clay M8-LC pas Fe 
‘ vat \ T Natural <a Ser oa Natural 
, eb ky 
i Na! \ 
4 g = 
> { Ws \ “Aba atti \ % 
2 200 4 | 4, g =, | @ 
i \ I ao Nevin at 
| 2 “a/ a 8 i oh 
\\ : a j a § i As vel etl : A 1 a y I 
ih’ ge amattont itil ‘ Bae 
ek: cn adstwdihentdil 
<6 pull au he i al 
iy Ah Aether 5 7 Be 28a VW a 
Neha Li pi i i“ Md ak eed hoch masochist poche Ved Nad 
T T T T T > av eB al Hd Tes tT : T T T T T T 1 
290 26 


Figure 7. Oriented clay diffractograms of the La Cantera Formation (Aptian) mudstone samples in stratigraphical order under three different 
conditions: calcined, glycolated, and natural. Main peaks identified are: I, illite; I/S, illite/smectite; K, kaolinite; S, smectite; and C, chlorite. 
1, Sample M1. 2, Sample M2. 3, Sample M3. 4, Sample M4. 5, Sample M5. 6, Sample M6. 7, Sample M7. 8, Sample M8. 
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under warm and humid conditions (Chamley, 1989). Chlo- 
rite characterizes all the samples as well; it is possible that 
along with illite, it is the result of smectite illitization and 
chloritization (Chamley, 1989). However, these processes 
characterize very high diagenetic stages, which is unlikely 
for our samples given the smectite presence in our samples. 
Therefore, the studied section was not submitted to signifi- 
cant deep burial and chlorite and illite are likely detrital in 
origin (Tineo et a/,, 2022). It is interesting that calcite is an 
abundant mineral in the lacustrine deposits, because it in- 
dicates that lake waters were alkaline. It is possible that the 
abundant calcite content is linked to the gypsum found in 
more marginal facies that were subaerially exposed (Castillo 
Elias, 2016). An alkaline paleolake is also confirmed by the 
presence of analcime in the clay mineral assemblage, since it 
precipitates in alkaline waters in association with smectite; 
commonly, those conditions occur when there is a source of 
detrital volcanics (Chamley, 1989; Salduondo et a/., 2020). 
Moreover, recent mineralogical results from Triassic paleo- 
lakes have also reported analcime related to high input of 
volcaniclastic material and alkaline waters (Benavente et al, 
2015; Salduondo et al, 2020). 


DISCUSSION 
Diagenetic overprint 

In order to accurately interpret the different stable iso- 
tope datasets obtained from the fossil remains tested, it is 
necessary to conduct a diagenetic assessment. Therefore, 
the correlation-regression analysis of 5’°O¢93 and 6'8Opq, 
was applied to fish remains. Since both components are 
formed from the same oxygen source, a positive correlation 
is expected between both 6'80 (Kolodny et al, 1983; lacumin 
et al, 1996; Tutken et al, 2006). In the case of the La Cantera 
fossil fish, the correlation coefficient found was r = -0.28, 
meaning that there is a negative correlation between 
5180-93 and 8180p, of 28%. This value is very low in com- 
parison with other values for published datasets, like r = 0.76 
obtained by Tutken et a/. (2006) or r = 0.98 from lacumin et 
al. (1996). Therefore, it is probable that the La Cantera fish 
remains analyzed present an altered isotopic signature to 
some extent. However, this does not assure that all values 
obtained have been diagenetically altered. As has been ex- 


tensively discussed in the literature, PO, exchange is negli- 
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gible in inorganic solutions (Kolodny et a/., 1983; lacumin et 
al, 1996). 

The sampled tissues, such as dermal bones and ele- 
ments associated with fins like fulcra, are covered by 
ganoine. Despite some authors proposing ganoine as a ho- 
mologue tissue of tooth enamel of tetrapods (Qu et al, 
2015), the most accepted idea is that ganoine is an exclusive 
feature of Actinopterygii (Schultze, 1996, 2016), being dif- 
ferent from true enamel considering its growth and histo- 
logical arrangement (Schultze, 2018). Therefore, different 
formulas were applied to our data to rule out biases due to 
the nature of the sampled materials. 


Paleotemperature calculations 

Regarding the paleotemperature calculations obtained, 
the values resemble what has been found in the diagenesis 
assessment, since most of the values are reasonable for a 
water column in a shallow paleolake in a general warm 
paleoclimate context for the Early Cretaceous (Ufnar et al, 
2004). Five of the values obtained are unreasonably warm 
(40-45 °C) for most fishes to inhabit, implying that those 
remains suffered diagenetical reset. These samples repre- 
sent the 30% of the data set and they don’t imply the 
method used is not reliable, but probably there are external 
factors for the preservation of the isotopic signature that 
require further investigation. 

Extant analogous ganoid Actinopterygii live up to 30 °C— 
e.g, the gar Lepisosteous oculatus (Grande, 2010)—. More- 
over, according to breeding places data, the bichires like 
Polypterus senegalus and P. de/hezi tolerate a range of tem- 
perature between 24-28 °C. Besides, species of the basal 
actinopterygians Order Acipenseriformes can tolerate up to 
20 °C, according to the Fish Base (Froese & Pauly, 2023). 
Alternatively, it is possible that during extremely low lake 
levels due to the dry season, temperatures became lethal 
for part of the biota. This has been documented in Lake 
Natron when lake levels descend to two or three m and 
water column temperatures rise to 35-45 °C. Considering 
other factors affecting the isotopic signature that might 
have led to paleotemperature anomalies (like physiology), 
the correlation-regression analysis was applied to each 
taxonomic group exclusively and to each tissue sampled, 


finding no difference in the results obtained. Nevertheless, 
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it is possible that the nature of the sampled material af- 
fected the preservation of the original isotopic composition 
and it has not been detected due to sample size. This might 
be a feasible possibility, since in Table 1 it is shown that the 
unreasonable paleotemperature values are from bones and 
some of the teeth, but particularly the ones that were at- 
tached to jaw fragments (bone). 

That 70% of our dataset does reflect reasonable pa- 
leotemperature calculations (23.33—35.80 °C) is consid- 
ered for the discussion as following. Within those values, 
there is no clear representation of a cold end member in the 
paleothermal range (minimum 23.33 °C). However, this 
would be expected in a warmer paleothermal context as has 
been described for the La Cantera Formation (Pramparo et 
al, 2018). Even more, the La Cantera underfilled system 
most likely was of restricted extent, with low thermal mass 
and prone to heat easily in an arid context, similar to what 
has been observed for modern lakes (Woolway & Merchant, 
2017). The paleotemperature values found are very similar 
to temperatures measured during the last decades for 
modern lakes in Africa and Asia in extensional tectonic 
context, like Lake Kivu (surface: 2,370 km2, T = 23 °C); Lake 
Turkana (surface: 7,785 km2, T = 28.5 °C); Lake Chilka (sur- 
face: 885 km?, T = 25 °C); and Lake Upemba (surface: 573 
km?, T = 26 °C) (Liu et a/, 2019). Paleotemperatures inter- 
preted as reasonable for the La Cantera water column are 
also in agreement with values found for other Cretaceous 
lacustrine deposits from USA, in which temperatures reach- 
ing 37.8 °C have been obtained from lacustrine carbonates 
(Fetrow et al, 2022). 

This implies that the stable isotope signature of fish 
remains could be used to some extent as a paleoclimate 
proxy coupled with other reliable paleoclimate proxies and 
applying diagenesis assessment to the data retrieved. Con- 
sistently with that, 6"3C values obtained from plant remains 
from the La Cantera Formation (average: -25%o) are re- 
markably similar to the average 6"3C of modern C3 plants 
(-25.4%o) (Craig, 1954), wood in particular (Grocke, 2002) 
which is consistent with the sampled stems. These strongly 
suggests that plant remains preserve the isotopic primary 
signature. However, the number of specimens that pass the 
preservation analysis was low, so preservation is a signifi- 


cant and limiting factor to consider to further test this proxy. 


Our dataset has an average approximate 5'C-97 composi- 
tion of -5.0%o &'°C. This matches previous values obtained 
for the Cretaceous of -3.0%o "°C (Grocke, 2002), supporting 
that with careful and extensive sampling, the La Cantera 
plant fossils would allow exploring regional climate variables 
such as 8'3C of atmospheric CO (Arens & Jahren, 2000; Grocke, 
2002). Furthermore, according to modern experiments on 
C3 plants, the values we have found would correspond to 
CO, enrichment conditions (Beerling & Woodward, 1995). 


Mineralogical proxies 

Regarding the clay mineral assemblage of the La Cantera 
mudstones, most of the analyzed samples are enriched in 
mixed-layer illite/smectite (I/S). This could be related with 
the progressive transformation of smectite to illite during 
burial (Hower et a/,, 1976). Therefore, it is possible that I/S 
is related with the presence of former smectite and conse- 
quently reflecting relatively warm-subarid conditions (alter- 
nation of wet and dry seasons) in the hinterlands (Chamley, 
1989). Overall, most of the samples from the studied sec- 
tion are dominated by illite and chlorite, and to a lesser ex- 
tent by kaolinite and I/S. Illite and chlorite are usually 
associated with relatively low chemical weathering forming 
through physical weathering (eg.,, Chamley, 1989; Hillier et 
al, 1995; Tineo et al, 2022). 

These results point to an increased hydrolysis index, 
meaning significant seasonal rainfall (Singer, 1984; Chamley, 
1989). Similar interpretations have been obtained from clay 
assemblages of the Cretaceous of Salta Province, north of 
Argentina (Do Campo et a/., 2010). A seasonal context is in 
agreement with a shallow lake that probably lowered its 
water column during the dry season, fluctuating from warm 
to warmer temperatures, like it is indicated by the paleotem- 
perature dataset obtained. The clay mineral assemblage 
interpretation is also in agreement with previous paleocli- 
mate interpretations made from palynological studies that 
found Classopolis and Gnetales specimens, indicating aridity 
episodes during La Cantera deposition (Pramparo et al, 
2018). Even more, there are several studies recording 
warm paleotemperatures and aridity conditions for mid- 
Cretaceous globally (Fetrow et a/., 2022) and particularly for 
terrestrial ecosystems during the Aptian—Albian (Ludvigson 
et al, 2015). 
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Another significant sedimentological feature of the La 
Cantera deposits is that it contains abundant primary gypsum 
deposits in other sections of the unit throughout the finely 
laminated mudstones. These have led to the interpretation 
of an underfilled lake basin type for the La Cantera lacus- 
trine system (Castillo Elias, 2016) sensu Bohacs et al, (2000), 
in which evaporites characterize the sedimentary column. 
This provides a coherent general context for the informa- 
tion we have found in this study. 


CONCLUSIONS 

The paleontological dataset from the La Cantera For- 
mation (fish and macroflora) is a useful and accurate mean 
of developing paleoclimate proxies. Fish remains preserve 
the primary isotopic signature partially. Therefore, coupled 
with other paleoclimate proxies and with rigorous diagene- 
sis tests, it can potentially be applied for paleoclimate re- 
constructions for this system. Plant remains remarkably 
preserve the original isotopic signature in the range of 
modern C3 plants and allow constraining paleoclimate and 
paleoenvironmental conditions. Furthermore, they consti- 
tute a potential proxy for deep time local C biogeochemical 
cycle reconstruction and global 6'7C-gz composition, since 
our results show agreement with other localities atmospheric 
CO, carbon isotopic composition during the Cretaceous. 
Therefore, these remains should be further investigated. 

The isotopic signature from unaltered fish remains al- 
lowed obtaining a paleothermal range of 23.33-35.8 °C for 
the water column of the La Cantera paleolake. These values 
are warm with a high minimum temperature, in agreement 
with the Lower Cretaceous global paleoclimate context in 
which maximum heat gain due to latent heat flux during a 
hothouse occurred at ~50 °N, at an equivalent paleolatitu- 
dinal range of the La Cantera lacustrine system in the 
Southern Hemisphere. Our data also agrees locally with an 
underfilled lake basin development in an intracontinental 
rift basin. Clay mineral assemblage composition indicates a 
seasonal rainfall distribution that supports the development 
of an underfilled paleolake that experienced seasonal low 
lake levels, warming up its water column. This is in agree- 
ment with previous palynological and sedimentological 


interpretations of aridity for the La Cantera Formation. 
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